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Abstract 

In this paper, the Lie symmetry method is performed for the nonlinear 
evolution equation for description nonlinear waves in a viscoelastic tube, 
we will find one and two-dimensional optimal system of Lie subalgebras. 
Furthermore, preliminary classification of its group-invariant solutions are 
investigated. 
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1 Introduction 

The theory of Lie symmetry groups of differential equations was developed by 
Sophus Lie [10]. Such Lie groups are invertible point transformations of both 
the dependent and independent variables of the differential equations. The 
symmetry group methods provide an ultimate arsenal for analysis of differential 
equations and is of great importance to understand and to construct solutions of 
differential equations. Several applications of Lie groups in the theory of differ- 
ential equations were discussed in the literature, the most important ones are: 
reduction of order of ordinary differential equations, construction of invariant 
solutions, mapping solutions to other solutions and the detection of linearizing 
transformations (for many other applications of Lie symmetries see [11], [2] and 

In the present paper, we study the following fifth-order nonlinear evolution 
equation 

u t + auu x + bu x 3 + cu x 4 + du x s = eu x 2. (1) 

where a 7 b,c,d and e are positive constants. This equation was introduced recently 
by Kudryashov and Sinelshchikov [8] which is the generalization of the famous 
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Kawahara equation. By using the reductive perturbation method, they obtained 
the equation (1). Study of nonlinear wave processes in viscoelastic tube is the 
important problem such tubes similar to large arteries (see [4] , [9] and [3] ) . 

In this paper, by using the lie point symmetry method, we will investigate 
the equation (1) and looking the representation of the obtained symmetry group 
on its Lie algebra. We will find the preliminary classification of group-invariant 
solutions and then we can reduce the equation (1) to an ordinary differential 
equation. 

This work is organized as follows. In section 2 we recall some results needed 
to construct Lie point symmetries of a given system of differential equations. 
In section 3, we give the general form of a infinitesimal generator admitted 
by equation (I) and find transformed solutions. Section 4, is devoted to the 
construction of the group- invariant solutions and its classification which provides 
in each case reduced forms of equation (I). 



2 Method of Lie Symmetries 

In this section, we recall the general procedure for determining symmetries for 
any system of partial differential equations (see [11] and [2]). To begin, let us 
consider the general case of a nonlinear system of partial differential equations 
of order nth in p independent and q dependent variables is given as a system of 
equations 

A»(i,u w )=0, v = l, ■•-,1, (2) 

involving x = (x 1 , ■ ■ ■ , x p ), u = (u 1 , ■ ■ ■ , u q ) and the derivatives of u with respect 
to x up to n, where vS n ^ represents all the derivatives of u of all orders from 
to n. We consider a one-parameter Lie group of infinitesimal transformations 
acting on the independent and dependent variables of the system (2) 

x { = x { + sC (x,u) + 0(s 2 ), i = l---,p, 

(3) 

v? = u J + srf (x, u) + 0(s ), j = 1 • • • , q, 

where s is the parameter of the transformation and , rf are the infinitesimals of 
the transformations for the independent and dependent variables, respectively. 
The infinitesimal generator v associated with the above group of transformations 
can be written as 

p <? 
v^^u^+^^i,^. (4) 
i=i j=i 

A symmetry of a differential equation is a transformation which maps solutions 
of the equation to other solutions. The invariance of the system (2) under the 
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infinitesimal transformations leads to the invariance conditions (Theorem 2.36 
of [11]) 

Pr w v[A„(a;,iiW)]=0, v = l,---,l, whenever A„(x,ii ( ")) = 0, (5) 

where Pr^ is called the n th order prolongation of the infinitesimal generator 
given by 

a = l J 

where J = ■ ■ ■ ,jk), 1 < jk _• P, 1 _• & _• i and the sum is over all J's of 
order < #J < n. If # J = k, the coefficient ipj of <9 U » will only depend on 
fc-th and lower order derivatives of u, and 

i=l i=l 

where := du a /dx l and Wj^ := duj/dx 1 . 

One of the most important properties of these infinitesimal symmetries is 
that they form a Lie algebra under the usual Lie bracket. 

3 Lie symmetries of the equation (1) 

We consider the one parameter Lie group of infinitesimal transformations on 

(x 1 = x, x 2 = t, u 1 = u), 

x = x + s£(x,t,u) + 0(s 2 ), 

t = x + sr)(x,t,u) + 0(s 2 ), (8) 
u = x + S(p(x,t,u) + 0(s 2 ), 

where s is the group parameter and £ x = £, £ 2 = rj and ip 1 — ip are the in- 
finitesimals of the transformations for the independent and dependent variables, 
respectively. The associated vector field is of the form: 

v = £(x, t, u)d x + rj(x, t, u)d t + <p(x, t, u)d u . (9) 

and, by (6) its third prolongation is 

Pr (5) v = v + <p x d Ux + <p* d Ut + <p?' d Ux2 + <p xt d Uxt + ■ ■ ■ 

• •• + ^X 2 +^* 4 a„ xt4 +/a„ t5 . (io) 
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where, for instance by (7) we have 

tp x = D x (tp-£u x -r)u t )+t;u x 2+riu xt , 
ip* = D t ((p~^u x ~r]Ut)+^u x t + r]u t 2 1 

(11) 

<P* = Dl(ip-^U x -7]U t )+^U x 5 t + 7]U t 5, 

where D x and D t are the total derivatives with respect to x and t respectively. 
By (5) the vector field v generates a one parameter symmetry group of equation 
(1) if and only if 

Pr^ 5 V[u 4 + auu x + bu x s + cu x a + du x & — eu x 2] = 0, 

(12) 

whenever ut + auu x + bu x 3 + cu x 4 + du x & — eu x 2 = 0. 
The condition (12) is equivalent to 

au x ip + auip x +ip t — e(p x2 + b(p x * + ap** + dip x& = 0, 

(13) 

whenever u t + auu x + bu x 3 + cu x 4 + du x a — eu x 2 = 0. 

Substituting (11) into (13), and equating the coefficients of the various mono- 
mials in partial derivatives with respect to x and various power of u, we can find 
the determining equations for the symmetry group of the equation (1). Solving 
this equation, we get the following forms of the coefficient functions 

£,=c 2 at + c 3 , T] = ci, ip = c 2 . (14) 

where c\, c 2 and C3 are arbitrary constant. Thus, the Lie algebra of infinitesimal 
symmetry of the equation (1) is spanned bye the three vector fields: 

V! = d x , v 2 = d t , v 3 =td x + ^d u . (15) 
The commutation relations between these vector fields are given in the Table 1. 



Table 1: The commutator table 



[Vi.Vj] 


Vl 


v 2 


V3 


Vl 











v 2 








Vl 


V3 





-vi 
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Theorem 3.1 The Lie algebra £3 spanned by vi,V2,v 3 is second Bianchi class 
type and it's solvable and Nilpotent. [7] 

To obtain the group transformation which is generated by the infinitesimal 
generators Vj for i = 1,2,3 we need to solve the three systems of first order 
ordinary differential equations 

dx(s) 



ds 
dt(s) 

ds 
du(s) 



= &(x(s),t(s),u(s)), x(0) = x, 
= Vl (x(s),t(s),ii(s)), t(0)=t, i = l,2,3 (16) 
= (pi(x(s),t(s),u(s)), u(0) = u. 



ds 

Exponentiating the infinitesimal symmetries of (1), we get the one-parameter 
groups Gi(s) generated by Vj for i = 1, 2, 3 



G 2 
G 3 

Consequently, 



(i, x,u) 1 — > (x + s,t,u), 

(t,x,u) 1 — > (x,t + s,u), (17) 
(t,x,u) 1 — > (x + ts,t,u + s/a). 



Theorem 3.2 If u = f(x,t) is a solution of (1), so are the functions 
Gi(«)-/(a:,t) = f(x-s,t), 

G 2 (s)-f(x,t) = f(x,t-s), (18) 
G 3 (s)-f(x,t) = f(x-ts,t) + s/a. 



4 Optimal system and invariant solution of (1) 

In this section, we obtain the optimal system and reduced forms of the equation 
(1) by using symmetry group properties obtained in previous section. Since the 
original partial differential equation has two independent variables, then this 
partial differential equation transform into the ordinary differential equation 
after reduction. 

Definition 4.1 Let G be a Lie group with Lie algebra q. An optimal sys- 
tem of s— parameter subgroups is a list of conjugacy inequivalent s— parameter 
subalgebras with the property that any other subgroup is conjugate to pre- 
cisely one subgroup in the list. Similarly, a list of s— parameter subalgebras 
forms an optimal system if every s— parameter subalgebra of q is equivalent to 
a unique member of the list under some element of the adjoint representation: 
F=Ad( 5 0>)).[ll] 
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Theorem 4.2 Let H and H be connected s- dimensional Lie subgroups of the 
Lie group G with corresponding Lie subalgebras f) and f) of the Lie algebra g of 
G. Then H^gHg^ 1 are conjugate subgroups if and only if f) = Ad(g(f))) are 
conjugate subalgebras. [11] 

By theorem (4.2), the problem of finding an optimal system of subgroups 
is equivalent to that of finding an optimal system of subalgebras. For onc- 
dimcnsional subalgebras, this classification problem is essentially the same as 
the problem of classifying the orbits of the adjoint representation, since each one- 
dimensional subalgebra is determined by nonzero vector in g. This problem is 
attacked by the naive approach of taking a general element V in g and subjecting 
it to various adjoint transformation so as to "simplify" it as much as possible. 
Thus we will deal with th construction of the optimal system of subalgebras of 
0- 

To compute the adjoint representation, we use the Lie series 

e 2 

Ad(exp(£v l )v i ) = Vj - e[vi,Vj] + —fa, [vi,vj]] , (19) 

where [vj,Vj] is the commutator for the Lie algebra, e is a parameter, and 
i,j = 1, 2, 3. Then we have the Table 2. 

Table 2: Adjoint representation table of the infinitesimal generators v» 



Ad(exp(evj))vj 


Vl 


v 2 


v 3 


Vl 


Vl 


v 2 


V3 


v 2 


Vl 


V 2 


v 3 - evi 


v 3 


Vl 


v 2 + evi 


v 3 



Theorem 4.3 An optimal system of one- dimensional Lie algebras of the equa- 
tion (1) is provided by 1) V2, 2) v 3 + av2 

Proof: Consider the symmetry algebra g of the equation (1) whose adjoint 
representation was determined in table 2 and 

V = oivi + a 2 v 2 + a 3 v 3 . (20) 

is a nonzero vector field in g. We will simplify as many of the coefficients a» 
as possible through judicious applications of adjoint maps to V. Suppose first 
that a 3 7^ 0. Scaling V if necessary we can assume that a 3 = 1. Referring to 
table 2, if we act on such a V by j4c?(exp(aiV2)), we can make the coefficient of 
vi vanish and the vector field V takes the form 

V = Ad(exp(aiv 2 ))V = a' 2 v 2 + v 3 . (21) 
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for certain scalar a' \. So, depending on the sign of a' 2, we can make the coeffi- 
cient of v 2 either +1, -1 or 0. In other words, every one-dimensional subalgebra 
generated by a V with 03 ^ is equivalent to one spanned by either V3 + V2, 
V3 - v 2 or v 3 . 

The remaining one-dimensional subalgebras are spanned by vectors of the 
above form with a 3 = 0. If a 2 ^ 0, we scale to make a 2 = 1, and then the vector 
field V takes the form 



for certain scalar a"\. Similarly we can vanish a"\ , so every one-dimensional sub- 
algebra generated by a V with 03 = is equivalent to the subalgebra spanned by 



Theorem 4.4 An optimal system of two-dimensional Lie algebras of the 
equation (1) is provided by 

< av 2 + v 3 , jSvi + 7V3 > 

Symmetry group method will be applied to the (1) to be connected directly to 
some order differential equations. To do this, a particular linear combinations of 
infinitesimals are considered and their corresponding invariants are determined. 

The equation (1) is expressed in the coordinates (x,t,u), so to reduce this 
equation is to search for its form in specific coordinates. Those coordinates 
will be constructed by searching for independent invariants (x, C) corresponding 
to the infinitesimal generator. So using the chain rule, the expression of the 
equation in the new coordinate allows us to the reduced equation. 

In what follows, we begin the reduction process of equation (1). 

4.5 Galilean-Invariant Solutions. First, consider V3 = td x + ^<9 M . To 
determine independent invariants /, we need to solve the first partial differential 
equations Vi(I)=0, that is invariants £ and x can be found by integrating the 
corresponding characteristic system, which is 



V" = 



a 1V1 + v 2 . 



(22) 



□ 



dt 
"0 



dx 



t 



a du 



1 



(23) 



The obtained solution are given by 



c 



x 



(24) 



X 



t 



Therefore, a solution of our equation in 



this case is 



u = f{x,x,0 




(25) 
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The derivatives of u are given in terms of £ and % as 

1 1 , N 

u x = — , u x 2 = u x 3 = u x 4 = u x 5 = 0, u t = C x - a;- (26) 

Substituting (26) into the equation (1), we obtain the order ordinary differential 
equation 

Cx + I C = 0. (27) 

A 

The solution of this equation is C = Consequently, we obtain that 

. . s + aci 
u(a?,t)= — - ■ (28) 

4.6 Travelling wave solutions. The invariants of V2 + cq vi = Co + 9t are 

X = x — c 1 and ^ = u so the reduced form of equation (1) is 

-coCx + a CC x + b( x z +c( x * +d( x 5 -eC x 2 = 0. (29) 

The family of the periodic solution for Eq.(29) when a — 1 takes the following 
form (see [8]) 

( = a a + A sn 4 {m x,k} + B sn{m \, k} -^-sn{m \, k}. (30) 

aX 

where sn{m x, k} is Jacobi elliptic function. 

4.7 The invariants of v 3 +/3v 2 = td x + f3d t + \d u arc \ = x ~ an d C — u ~ ^3 
so the reduced form of equation (1) is 



1 t 



— - -Cx + «CCx + K x 3 + C C X 4 +d( x 5 - e (x 2 = °- ( 31 ) 



4.8 The invariants of v 2 = d t arc x — x an d C = u then the reduced form of 
equation (1) is 

aCCx + Kx 3 + C C X 4 +^Cx s ~ e C x 2 = °- ( 32 ) 

4.9 The invariants of vi = d x arc x = t and ( = u then the reduced form of 
equation (1) is ( x = 0, then the solution of this equation is u(x, t) = cte. 

Acknowledgment 

This research was supported by Islamic Azad University of Karaj braench. 



8 



References 



[I] G.W. Bluman and J.D. Cole, Similarity Methods for Differential Equa- 
tions, Applied Mathematical Sciences, No. 13, Springer, New York, 1974. 

[2] G.W. Bluman and S. Kumei, Symmetries and Differential Equations, 
Springer, New York, 1989. 

[3] C.G. Caro, T.J. Pedly, R.C. Schroter and W.A. Seed, The mechan- 
ics of the circulation, Oxford: Oxford University Press, 1978. 

[4] Y.C. Fung, Biomechanics: mechanical properties of living tissues, New 
York, Springer- Verlag, 1993. 

[5] I.L. Freire and A.C. Gilli Martins, Symmetry coefficients of semilinear 
PDEs, arXiv: 0803.0865vl, 2008. 

[6] N.H. IBRAGIMOV, (Editor), CRC Handbook of Lie Group Analysis of Dif- 
ferential Equations, Vol. 1, Symmetries, Exact Solutions and Conservation 
Laws, CRC Press, Boca Raton, 1994. 

[7] S. V. Khabirov, Classification of Three- Dimensional Lie Algebras in R 3 
and Their Second- Order Differential Invariants, Lobachevskii Journal of 
Mathematics, 31 (2) (2010), 152-156. 

[8] N.A. KUDRYASHOV and D.I. Sinelshchikov, Nonlinear evolution equation 
for describing waves in a viscoelastic tube, Commun Nonlinear Sci Numer 
Simulat, 16 (2011), 2390-2396. 

[9] N.A. KUDRYASHOV and I.L. Chernyavskii, Nonlinear waves in fluid flow 
through a viscoelastic tube, Fluid Dynam, 41(1) (2006), 49-62. 

[10] S. Lie, Theories der Tranformationgruppen, Dritter und Letzter Abschnitt, 
Teubner, Leipzig, 1893. 

[II] P.J. Olver, Applications of Lie Groups to Differential Equations, 
Springer, New York, 1993. 

[12] A.D. POLYANIN, and V.F. ZAITSEV, Handbook of Nonlinear Partial Dif- 
ferential Equations, Chapman & Hall/CRC, Boca Raton, 2004. 



9 



